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INTRODUCTION 

The helicopter rotor design process is highly multidisciplinary in nature 
and requires a merging of several technical disciplines such as dynamics, 
aerodynamics, structures and acoustics. In the past the conventional design 
process was controlled by the designer's experience and the use of trial and 
error methods. Today, one of the more promising approaches to the rotor 
blade design process is the application of structural optimization 
techniques. An extensive amount of work has been done in developing design 
optimization procedures to bring the state of the art to a very high 

level 1 *^. While these techniques have received wide attention in the fixed- 

wing field 1 , they are fairly recent in the rotary wing industry^ -5 . Most of 
the work involving application of optimization techniques to rotor blade 
design has been focused on nearly independent technical disciplines with 
very little consideration of the coupling and interaction between the 
disciplines. For example, the dynamic design requirements have been 
considered in the optimum rotor blade design in refs. 6-10. Blade 
aerodynamic and structural requirements were considered in refs. 11 and 12, 
respectively . 

The necessity of merging appropriate disciplines to obtain an integrated 
design procedure has been recently emerging and with improved understanding 
of helicopter analyses and optimization schemes, it is now possible to apply 
optimization techniques and include the couplings between the disciplines. 
In refs. 13-15 the dynamic and structural design requirements were coupled 
with airloads in the analysis and in refs. 16 and 17 the dynamic and 
aeroelastic requirements were integrated. The optimization procedure 

described in this paper is part of an effort at NASA Langley Research 
1 8 

Center and is aimed at integrating two technical disciplines, aerodynamics 
and dynamics. As a first investigation, the airloads will be included to 
perform coupled airload/dynamic integration of rotor blades. Later the 
aerodynamic performance requirements will be added to obtain an integrated 
aerodynamic /dynamic optimum design procedure. The procedure is no longer 
sequential - rather it will account for the interactions between the two 
disciplines simultaneously. The paper briefly describes some of the recent 
work done by the authors which focussed on optimum blade design with dynamic 
behavioral constraints and presents some of the authors' recent experiences 
in developing a strategy for structural optimization with integrated 
dynamics/aerodynamics of rotor blades. 
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INTEGRATED ROTORCRAFT ANALYSIS 

Currently at the NASA Langley Research Center, there is an effort to 
integrate various technical disciplines such as dynamics, aerodynamics and 
structures into the rotor design process. Shown below in fig. 1 is a 
tentative plan of the integrated rotor analysis program. The plans are to 
perform independent discipline level optimizations, (e.g. rotor aerodynamic, 
dynamic and structural optimization as shown by the clear bubbles) by 
considering design variables, constraints and objective functions that 
affect the particular discipline considered. The next step is to couple 
rotor aerodynamics and dynamics to perform integrated aerodynamic/dynamic 
optimization. This would involve considerations of design variables and 
requirements of importance to each discipline, although there are certain 
design variables that influence all the disciplines involved. The 
structural design criteria are then introduced to obtain an integrated 
aerodynamic/dynamic/structural optimization procedure. The influence of 
airframe dynamics and acoustics will be accounted for through constraints in 
the design optimization to obtain the 'fully integrated procedure.' The 
final step is to validate this optimization procedure for a blade test 
article . 
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ROTOR BLADE DESIGN CONSIDERATIONS 

Rotor blade design involves several considerations some of which are listed 
below in fig. 2. The blade design must satisfy specified strength criteria 
and should be damage tolerant. The rotor blade aerodynamic design process 
consists of proper selection of blade geometric variables such as planform, 

l 1 

airfoils, twist, etc. to meet performance requirements ~. Helicopter 
performance is usually expressed in terms of horsepower required as a 
function of velocity. The horsepower required to drive the main rotor for 
any part of a mission must be less than the available horsepower. The 
airfoil section stall must also be avoided, i.e. the airfoil sections must 
operate at section drag coefficients less than a specified value. Two other 
major criteria in rotor blade design have been low weight and low vibration. 
For a helicopter in forward flight, the nonuniform flow passing through the 
rotor causes oscillating airloads on the rotor blades . These loads in turn 
are translated into vibratory shear forces and bending moments at the hub. 
Therefore, vibration alleviation without weight penalty is an important 

17 19 

criterion. The blade should also be aeroelastically stable ' and finally 
the noise levels generated by the rotor which are a function of local Mach 
number and airloads should be reduced. This paper will concentrate on the 
low vibration and the low blade weight aspects of the design. 


• Strength, survivability, fatigue life 

• Aerodynamic performance 

• Vibration 

• Weight 

• Aeroelastic stability 

• Acoustics 


FIGURE 2 
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DYNAMIC OPTIMIZATION PROBLEM STATEMENT 


As mentioned before, low vibration is an important design requirement in 
helicopter rotor blade design. One way of reducing the vibration level in 
the blade is to design the blade such that the natural frequencies are 
separated from multiples of the driving frequencies. Failure to consider 
frequency placement early in the design process can cause a significant 
increase in the final blade weight later if postdesign addition of 
nonstructural _ masses is required. Appropriately placing the natural 
frequencies can be done by a proper tailoring of the blade mass and/or 
stiffness distributions to meet the necessary design requirements using 
structural optimization. This section of the paper presents an overview of 
the dynamic optimization work which has been completed. The goal of the 
dynamic optimization problem (fig. 3) is to obtain minimum weight designs of 
blades with constraints on multiple coupled flap-lag natural frequencies. 
It is also important that the autorotational performance of the blade not be 
degraded during the tailoring process since the blade should have sufficient 
inertia to autorotate in case of an engine failure. In order to ensure a 
safe design, the blade centrifugal stress should be limited by an 
appropriate upper bound. For this study only centrifugal stress has been 
considered. The blade is assumed to be in vacuum in this investigation and 
the results of this analysis will generate a good starting point for the 
integrated optimization. 


• Goal - Minimize blade weight with constraints 

on multiple coupled natural frequencies, 
autorotational inertia and stress 

• Approach - Stiffness and/or mass modifications, 

placement of tuning masses 

• Assumption - Blade is in vacuum - generates 

good starting point for integrated 
optimization 


FIGURE 3 
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ROTOR BLADE MODEL FOR DYNAMIC OPTIMIZATION 


The rotor blade model for dynamic optimization is shown below in fig. 4. 
The blade is articulated and has a fixed hub, a pretwist and a root spring 
which allows torsional motion. A box beam with unequal vertical wall 
thicknesses is located inside the airfoil and lumped nonstructural masses 
are located inside the box and distributed spanwise. This model is based on 
an existing blade design denoted the 'reference blade' described in refs. 8, 
9, and 13. As in ref. 13, it is assumed that the box beam contributes all 
the blade stiffness, that is, the contributions of the skin, honeycomb, etc. 
to the blade flap and lag stiffnesses are neglected. The details for 
calculating the box beam section properties can be found in ref. 8. The 
properties of the box beam located inside the airfoil are as follows: 

h-0.117 ft, b-0.463 ft, p-8.645 slugs/ft 3 , E-2.304xl0 9 lb/ft 2 , allowable 

7 2 

stress a_- -1.93x10 lb/ft and factor of safety, FS-3. The blade is 

discretized into ten segments. Both rectangular and tapered blades are 
considered. For the rectangular blade, the box beam outer dimensions along 
the blade span remain unchanged. The design variables for the rectangular 
blade are the box beam wall thicknesses tj_, t 2 , and tj and the magnitudes of 

the nonstructural weights located inside the box beam at ten spanwise 
locations. For the tapered blade it is assumed, as in refs. 8 and 9 that 
the box beam is tapered and the additional design variables are the box beam 
height at the root, h r , and the taper ratio, which is defined as the 

ratio of the box beam height at the root to the corresponding value at the 
tip. A linear variation of the box beam height, h, in the spanwise 
direction is assumed. 


Lumped 
mass - 




= 


Reference blade k b 1 H h f /' "W j 

• Articulated, rigid hub JAs / 

• Rectangular planform, pretwist, / 

root spring b r \4 

Design variables ^ 

• Box beam wall thicknesses, \ v t 2 , t 3 (10 spanwise positions) 

• Box beam outer dimension h r 

• Taper ratio A,h . . . 

• Magnitudes of lumped masses (10 spanwise positions) 


>r^f 


FIGURE 4 
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The purpose of the optimization procedure, as described in fig. 5 below, is 
to minimize the weight W of the rotor blade while constraining the natural 
frequencies f^ to be within specified 'windows' (upper and lower bounds). 

An existing blade which is being used in a production helicopter has been 
selected as a baseline blade and will be referred to as the 'reference 
blade' . A modal analysis of the reference blade showed that the frequencies 
of interest were not near the n per rev (critical values) values where n 
denotes the total number of blades. Hence it was decided to define 
constraints to force the frequencies of the optimum blade to be close to 
those of the reference blade. The concept of 'windows' has been used since 
the nonlinear programming method used in this work cannot handle equality 
constraints. These windows, denoted by f v and f v (for the lower bound and 

*L *U 


upper bound on frequency, respectively) , are on the frequencies of the first 
three lead-lag dominated modes and the first two flapping dominated modes 
(elastic modes only) . The frequency windows are carefully selected to 
alleviate any shear amplification problem. A prescribed lower 


limit a on 


the blade autorctational inertia AI and an upper bound on the blade 


centrifugal stress 


c k have also been used. 


max 

Side constraints <J>.- 


and 


*i 


U 


(lower and upper bounds on the i th design variable <t> i ) have been imposed on 
the design variables to avoid impractical solutions. 


• Objective function 

• Minimum blade weight W 

W = W b + W 0 

• Constraints 

• Frequency windows on first 3 lead-lag 

and first 2 flapping elastic modes 

V , * 5 , i' U k=,A3A5 

• Lower bound on autorotational inertia 
Ai > a 


Upper bound on centrifugal stress 

a k * FS - %, ax 

Bounds on design variables 


<t> 


i 


* ♦ 


'u 


FIGURE 5 
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METHODOLOGY FOR DYNAMIC OPTIMIZATION PROBLEM 

The procedure described in this paper uses the program Comprehensive 

Analytical Model of Rotorcraft Aerodynamics and Dynamics (CAMRAD) 2< \ The 
modal analysis portion of the program CAMRAD which uses a modified Galerkin 
2 1 

approach has been used for the dynamic optimization problem. According to 
ref. 22, this approach is the preferred method for computing mode shapes and 
frequencies of structures having large radial variations in bending 

stiffness. The general purpose optimization program CONMIN 23 which uses the 
nonlinear programming method of feasible directions has been used for the 
optimization. The method of solution described below (fig. 6) starts with 
discretizing the blade into finite segments. In the search for the optimum 
vector of new design variables, CONMIN requires derivatives of the objective 
function and constraints. The user has the option of either allowing CONMIN 
to calculate derivatives by using forward differences, or by supplying those 
derivatives to CONMIN. In the work presented in this paper, the latter 
approach has been used. Analytical expressions for the derivatives of the 
objective function and the autorotational inertia constraint have been 
obtained. A central difference scheme has been used for the derivatives of 

the frequency constraints. The initial attempt 8,9 using forward differences 
gave highly inaccurate derivatives. 

The optimization process generally requires many evaluations of the 
objective function and the constraints before an optimum design is obtained. 
The process therefore can be very expensive if exact analyses are made for 
each evaluation. To reduce computational requirements, the optimization is 
based on the use of approximate analyses. In the present paper a piecewise 
linear analysis, based on first order Taylor Series expansions, is used. 
The approximate analyses should produce accurate characteristics of the real 
problem in a neighborhood of the current design which is continuously 
updated during optimization. The method has been found to be effective in 
the past (e.g., ref. 24) for providing accurate approximations. 

• Codes used 

•CAMRAD - Blade modal analysis (modified 

Galerkin approach) 

• CONMIN - Optimization (nonlinear programming 

approach - method of feasible directions) 

• Method of solution 

• Discretize the blade (10 finite segments) 

• Compute mode shapes and frequencies 

• Perform sensitivity analysis 

• Analytical derivatives of objective function, 

autorotational inertia constraint and 
stress constraints 

• Central differences for frequency constraint 

.derivatives 

• Use approximate analysis techniques 

FI6URE 6 
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DYNAMIC OPTIMIZATION RESULTS FOR RECTANGULAR AND TAPERED BLADES 

Results obtained by applying the dynamic optimization procedure to the 
design of both rectangular and tapered rotor blades are summarized here 
(fig. 7) . The table below depicts some of the representative results for 
the rectangular and tapered blades. For the rectangular blade the 40 design 
variables are the box beam wall thicknesses (t^, t 2 , t^) and the magnitudes 

of the nonstructural masses at ten spanwise locations. For the tapered 
blade with 42 design variables, the two additional design variables are the 
box beam height at the root and the taper ratio. In each table, column 1 
represents the reference blade data; column 2 gives the corresponding 
information for the optimum design for the rectangular blade with 
constraints on the five frequencies, autorotational inertia and stress; and 
column 3 gives results for the optimum design for the tapered blade with the 
same set of constraints. In all cases convergence to optimum designs 
typically has been achieved in 8-10 cycles. 

The table indicates that the optimum rectangular blade is 4.7 percent 
lighter than the reference blade and the optimum tapered blade is 6.2 
percent lighter than the reference blade. Although the first lead-lag 
frequency (f^) is at its prescribed upper bound after optimization, both 

frequencies are satisfactory as far as the shear amplification problem is 
concerned. The autorotational inertia constraint is also active (i.e, 
exactly satisfied) in all the cases. 




Optimum blade 


Reference 

blade 

Rectangular 
(40 design 
variables) 

Tapered 
(42 design 
variables) 

X h 
f 1 ,Hz 

1.0 

12.285 

1.0 

12.408* 

1.49 ! 

12.408* 

N 

X 

16.098 

16.075 

16.066 

f 3 , Hz 

20.913 

21.081 

20.888 

N 

X 

H- 

34.624 

34.823 

34.678 

f 5 * Hz 

35.861 

35.800 

35.507 

Autorotational ^ 
inertia(AI), lb-ft 2 

517.3* 

517.3* 

517.3* 

Blade weight, lb 

98.27 

93.61 

92.16 

Percent reduction 
in blade weight ! 

— 

4.74 

6.21 


! From reference blade • Active 


FIGURE 7 
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OPTIMUM HORIZONTAL WALL THICKNESS (t 1 ) DISTRIBUTIONS 
WITH MULTIPLE FREQUENCY AND STRESS CONSTRAINTS 

The optimum box beam horizontal wall thickness (t^) distributions along the 

blade span are shown below in fig. 8 and are compared with the corresponding 
distribution of the reference blade. On the left, the optimum distribution 
corresponds to the rectangular blade with 40 design variables (column 2, 
fig. 7) . On the right, the optimum distribution corresponds to the tapered 
blade with 42 design variables (column 3, fig. 7) . In both cases the 
optimum blade has a larger value of t^ than the reference blade at the blade 

tip. The explanation for this is as follows. The autorotational inertia 
can be increased with an increase in the moment arm and, therefore, the 
constraint on the autorotational inertia is satisfied easily if more mass is 
moved to the blade tip. However, the presence of the centrifugal stress 
constraint counteracts this tendency. Therefore, the net result is more 
blade mass towards the outboard region of the blade (although, not 
necessarily all at the tip) . 


— Reference 

— Optimum 



FIGURE 8 
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OPTIMUM VERTICAL WALL THICKNESS (tj) DISTRIBUTIONS 
WITH MULTIPLE FREQUENCY AND STRESS CONSTRAINTS 

The optimum box beam vertical wall thickness (t 2 > distributions along the 

blade span are shown below in fig. 9 and are compared with the corresponding 
distribution of the reference blade. On the left, the optimum distribution 
corresponds to the rectangular blade with 40 design variables (column 2, 
fig. 7) . On the right, the optimum distribution corresponds to the tapered 
blade with 42 design variables (column 3, fig. 7) . In both cases the 
optimum blade has a larger value of t 2 than the reference blade at the blade 

tip due to the presence of the autorotational inertial constraint as 
explained in the previous chart. However, the difference in magnitude 
between the optimum and reference blade value at the blade tip is not as 
significant as it is for the horizontal wall thickness t^. The nature of 

the horizontal and vertical wall thicknesses (t^ and t 2 / respectively) are 

also different as the former primarily affects the flapping frequency and 
the later affects the lead-lag frequency. 


— Reference 
— - Optimum 


.012 

.010 

Vertical wall - 008 
thickness (to) .006 

ft .004 

.002 


.000 


, 40 design variables 


"» 1 


( 


42 design variables 


i — i — — - 


.83 4.4 8.8 13.2 17.6 22. .83 4.4 8.8 13.2 17.6 22. 

Blade radius, ft Blade radius, ft 


Rectangular blade 


Tapered blade 


FIGURE 9 
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OPTIMUM NONSTRUCTURAL SEGMENT WEIGHT DISTRIBUTIONS 
WITH MULTIPLE FREQUENCY AND STRESS CONSTRAINTS 

Shown below (fig. 10) are the optimum and the reference blade nonstructural 
segment weight distributions along the blade radius for both the rectangular 
blade with 40 design variables (column 2, fig. 7) and the tapered blade with 
42 design variables (column 3, fig. 7) . For the rectangular blade (left 
side of the figure) the optimum blade has lower nonstructural weight 
throughout the blade span. However, for the tapered blade (right side of 
the figure) the optimum blade has larger nonstructural weight towards the 
blade tip than the reference blade. This is because the tapered blade has 
reduced structural weight requirements at the blade tip. Hence, in order to 
satisfy the autorotational inertia constraint, the nonstructural weight at 
the tip must increase. Even so the total weight of the optimum blade is 
still lower than that of the reference blade. 



Blade radius, ft 


— Reference 
— - Optimum 

42 design variables 




.83 4.4 8.8 13.2 17.6 22. 

Blade radius, ft 


Rectangular blade 


Tapered blade 
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STRATEGY AND TASKS FOR STRUCTURAL OPTIMIZATION 
WITH INTEGRATED DYNAMICS /AERODYNAMICS 

The structural optimization of helicopter rotor blades with integrated 
dynamics/aerodynamics involves both dynamic, aerodynamic and structural 
design variables, constraints and objective functions along with the blade 
dynamic/aerodynamic/structural analysis. Together with calculations of the 
associated sensitivity derivatives this can make the integrated optimization 
process very complicated and expensive. As a first step towards integrating 
dynamics and aerodynamics, it was decided to separate the aerodynamic 
effects into two parts: airloads and performance (fig. 11). The initial 

step in integrated dynamic /aerodynamic optimization will combine airloads 
and dynamics. The second step would involve addition of aerodynamic 
performance to obtain a fully integrated structural optimization procedure 
with dynamics/aerodynamics . The inclusion of airloads would allow 
calculation of hub shears and moments which enter into the objective 
function and/or constraints. This would allow the inclusion of blade 
aeroelasticity through either limits on the hub loads or the blade stability 
margin. The aerodynamic analysis would include trimming of the blade at 
each step of the design process for a specified flight condition. The trim 
analysis is in fact a coupled dynamic/aerodynamic/structural procedure. 

The integrated design process would require the use of more than one 
objective function in the design formulation. This is because it is 
difficult to single out an objective function as the primary requirement in 
an engineering system as complex as the rotor blade. This leads to the 
necessity of using multiple objective function techniques to formulate the 
optimization problem. Therefore, various multiple objective function 
techniques are- being investigated and a method called 'Global Criteria 
25 

Approach' is being examined. 


• Dynamic/aerodynamic/structural design 

variables and constraints 

• Include airloads first - integrated dynamic/airload 

optimization procedure 

• Add aerodynamic performance next - fully 

integrated dynamic/aerodynamic 
optimization procedure 

• Coupled trim analysis 

• Several objective functions - multiple 

objective function handling capability 
required . 

• Evaluate ’Global Criteria* approach for multiple 

objective optimization 


FIGURE 11 
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ANALYSIS COUPLINGS FOR STRUCTURAL OPTIMIZATION 
WITH INTEGRATED DYNAMICS /AIRLOADS 

Below is a schematic diagram that shows the general flow of information 
between the three major analyses involved in integrated airloads/dynamic 
optimization. Note that the three major disciplines are internally coupled. 
For instance, the blade aerodynamic analysis provides the airloads and 
control settings which are fed into the blade dynamic analysis. The blade 
dynamic analysis, based on this information, provides the blade natural 
frequencies, mode shapes, hub shears, moments, etc. If unsteady 
aerodynamics is included, the dynamic and aerodynamic analyses are coupled 
as shown by the dotted line in fig. 12 below. The information obtained from 
the dynamic analysis (shears/bending moments) are fed into the structural 
analysis box along with the airloads from the aerodynamic analysis to 
perform the trim analysis. The structural analysis is also used to compute 
the blade centrifugal stresses which are incorporated as constraints in the 
optimization process. 



Airloads 


FIGURE 12 
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COMPUTATIONAL CONSIDERATIONS FOR STRUCTURAL OPTIMIZATION 
WITH INTEGRATED DYNAMICS /AIRLOADS USING CAMRAD 

Some of the computational considerations involved in the structural 
optimization procedure with integrated dynamics /airloads is described below 

in fig. 13. The program CAMRAD is used for the aerodynamic and dynamic 
analyses of the rotor blade in forward flight . The program has been found 

o q o <r 

to be very reliable for analysis of helicopter rotors ' It uses a 

lifting line or blade element approach to calculate the section loading from 
the airfoil two-dimensional aerodynamic characteristics with corrections for 

22 

yawed and three-dimensional flow effects . The program also has the 
provision for including unsteady aerodynamics. 

Each intermediate design should satisfy the trim condition. The program 
CAMRAD offers two broad categories of trimming - the free flight case and 
the wind tunnel case. In the free flight case, the entire helicopter is 
trimmed to force and moment equilibrium whereas in the wind tunnel case the 
isolated rotor is trimmed to a prescribed operating condition. It is 
possible to use a free flight trim option for an isolated rotor in a wind 
tunnel since the trim option and the degrees for freedom representing the 
aircraft can be specified independently. However, the wind tunnel trimming 
options are more typical of a rotor in a wind tunnel without consideration 
of the complete rotorcraft. The wind tunnel trim option is selected for 
this analysis since the model used in this study is a wind tunnel model of a 
rotor. The trim option consists of trimming the rotor lift, drag and 
flapping angle with collective pitch, cyclic pitch and shaft angle. 


• Aerodynamic loads (forward flight) 

• Lifting line theory to calculate section 

loading from airfoil 2-D aerodynamic 
characteristics 

• Corrections for yawed and 3-D flow effects 

• Trim analysis 

• Wind tunnel trim for isolated rotor 

• Lift, drag and flapping angle with 
collective pitch, cyclic pitch and 
shaft angle 


FIGURE 13 
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FORMULATION OF THE STRUCTURAL OPTIMIZATION PROBLEM 
WITH INTEGRATED DYNAMICS /AIRLOADS 

The optimization problem addressed here uses blade weight and blade root 4 
per rev vertical shear as the objective functions to be minimized. The 
constraints are 'windows' on the coupled flap-lag natural frequencies to 
prevent them from falling into the critical ranges, a prescribed lower bound 
on the blade autorotational inertia and a maximum allowable upper bound on 
the blade stress. The design variables (fig. 14) are the blade spanwise 
stiffness distributions (El's and GJ) , the magnitudes of the lumped 
nonstructural masses distributed spanwise, the blade taper ratio and the 
root chord as shown below in the figure. The nonstructural masses which 
were used for frequency placement in the dynamics work discussed earlier 
will now be used for both frequency tuning as well as hub shear alleviation. 


Objective function: Blade weight and blade root 

vertical shear 

Frequencies, autorotational 
inertia, blade stress 

Stiffness and mass distributions, 
magnitudes of lumped/tuning 
masses, taper ratio, root chord 


Constraints: 
Design variables: 



't 


Jt= 


c r : Root chord 
Cf. Tip chord 
X : Taper ratio 


FIGURE 14 
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FLOWCHART OF THE STRUCTURAL OPTIMIZATION PROCEDURE 
WITH INTEGRATED DYNAMICS/ AIRLOADS 

The optimization procedure shown in the flowchart below (fig. 15) is 
initiated by identifying the blade preassigned parameters which are the 
parameters that are held fixed during optimization. The next step is to 
initialize the design variables and perform the internally coupled blade 
analysis which comprises blade aerodynamic, dynamic and structural analyses. 
A sensitivity analysis is part of the procedure and consists of evaluations 
of the derivatives of the objective function and the constraints with 
respect to the independent design variables. Once the sensitivity analysis 
is performed the approximate model is defined based on a standard 
approximation technique. Using CONMIN along with the approximate model 
updated design variable values are obtained. The process continues until 
convergence is achieved. 



FIGURE 15 
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SENSITIVITY ANALYSIS 


The conventional approach for performing sensitivity analysis is to 
calculate the derivatives either analytically or by using finite difference 
schemes. Since analytical expressions are seldom available and use of 
finite difference schemes is usually expensive and sometimes inaccurate, a 
9 7 

new method for obtaining the system sensitivity has been considered for 
the present work. The method enables one to calculate the sensitivity 
derivatives of the system solution with respect to a design variable from a 
set of simultaneous equations which are known as Global Sensitivity 
Equations (GSE) . In fig. 16 the system sensitivity equations are described 
in terms of a coupled system consisting of the boxes A, D, and S 
representing aerodynamics, dynamics and structures. Each discipline box is 
regarded as a set of mathematical operations that solves one of the sets of 
governing equations on the right to produce an output denoted by Y. For 
example, Y A denotes the output of the aerodynamic analysis. The coupling of 

the system is demonstrated in the figure below. The design variables are 
denoted by X. The quantities X and Y are in general vectors. Furthermore 
the subset of Y A entering D may be different from the subset of Y A entering 

S, although the subsets may overlap. 

Using chain rule on the governing equations as in ref. 27, the system 
sensitivity equations are derived. The sensitivity derivatives appear as 
the vector of unknowns. The coefficient matrix consists of partial 
derivatives of the output of the various disciplinary responses with respect 
to each other positioned off the diagonal and identity submatrices along the 
diagonal. Nonzero values of these partial derivatives reflect system 
couplings. The right hand side vector contains the partial derivatives of 
the disciplinary outputs with respect to a particular design variable (e.g. 
X^) . The coefficient matrix needs only to be formed and factored once for a 

given system and then back substituted using a new right hand side vector 
for every new design variable. Thus the method enables the computations of 
derivatives of complex internally coupled systems without having to perform 
expensive finite difference derivatives based on the entire system analysis. 



Coupled rotor blade analysis 


Global sensitivity equations: 


Governing equations: 

A ((X, Y d , Y S ), Y a ) = 0 
D((X > Y A> Y S ),Y D ) = 0 
S((X, Y d> Y a ), Ys) = 0 


Sensitivity derivatives 
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FIGURE 16 
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MULTIPLE OBJECTIVE FUNCTION - GLOBAL CRITERIA APPROACH 


As indicated before, the current optimization procedure requires a multiple 
objective function approach. Several methods have been proposed for the 
solution of multiobjective optimization problems. However, many of these 
methods suffer from a need for assigning relative priorities to the 
individual objective functions, e.g. assigning weight factors. The 
optimization goal is to find the set of design variables <t> which minimizes N 
objective functions (F^ (<()), F 2 ($), . .., F N (<(>)) subject to a set of 

inequality constraints g^ ( j**l , 2 , . . . ,NCON where NCON denotes the total 

number of constraints) . Using the Global Criteria Approach described in 

★ 

fig. 17, the optimum solution $ is obtained by minimizing a prescribed 
'global criterion' f ($) which is defined as the sum of the squares of the 
relative deviations of the individual objective functions F^ (<|>) from their 

★ ★ 

respective feasible optimum values F^(<J>^). The optimum solution, <|>^, to the 

i t ^ 1 individual objective function is obtained by minimizing F^($) subject to 

the constraints g . (<(>) <0 , j-1, 2, . . . , NCON . The optimization problem now is to 

J ./•» 

minimize the composite objective function F(<j>) subject to exactly the same 
set of constraints as used in the individual optimizations. The method is 
less judgmental in the sense it imposes equal priority to each individual 
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objective function . 


Optimization goal 

Minimize "N" objective functions 

F 1 (♦)> F 2 (♦)■ F 3 (<►)>•••» f n (♦) 
subject to gj (<|>) <0 j = 1,2,..., NCON 

Global criterion formulation 

Minimize r_ ... - .. i2 


N 

F W=i!l 


F j W “ F j(4> j ) 




F|(+|) 

subject to g . (<|>) < 0 j = 1 ,2,..., NCON 

($]) obtained from 
Minimize Fj (<|>) 

subject to 9 j (<t>) < 0 j = 1 , 2 ,..., NCON 


FIGURE 17 
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FORMULATION OF STRUCTURAL OPTIMIZATION PROBLEM 
WITH INTEGRATED DYNAMICS /AIRLOADS USING GLOBAL CRITERIA APPROACH 

Using the the Global Criteria Approach the airload/dynamic optimization 
problem with multiple objective functions can be formulated as shown in fig. 
18. The two objective functions F-^ (<|> ) and F 2 (<t>) are the blade weight W and 

the blade root 4 per rev vertical shear F z , respectively. The constraints 

are on the frequencies f^, k**l,2,...6 (three lead- lag and three flapping 

dominated modes) , the blade stress a and the blade autorotational inertia 
AI. Using the Global formulation the new global objective function F(<|>) is 

defined as the sum of the squares of the deviations of the objective 

★ 

functions, W and F z , from their respective individual optimum values W and 

f£. The optimization problem now is to minimize F(<(>) subject to the 
original set of constraints. 


Multiple objective functions: F 1 (4>) = W 

F 2 (4>) = F z 


Constraints, g(<|>): 


Global objective function: 

a f\N - W 

m = v - 

subject to gfijO'SO 


1 • f k/ f|< L - 0 

f k/f ku -i 5 0 

a “ Al ^ 0 
a- FS- o max - 0 



FIGURE 18 
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STUDY OF GLOBAL CRITERIA APPROACH FOR HEIGHT-STRESS OPTIMIZATION 

(BLADE IN VACUUM) 

Before attempting to solve the above integrated airload/dynamic optimization 
problem it was first decided to study the Global Criteria Approach for the 
dynamic optimization problem with the blade in vacuum and the blade weight 
and centrifugal stress as the two objective functions to be minimized (fig. 
19) . There F^ is equal to W which is the blade weight and F2 is equal to a 

which represents the maximum centrifugal stress in the blade. The 
constraints are windows on the first coupled lead-lag dominated and the 
first flapping dominated frequencies and the blade autorotational inertia. 
The formulation of the test problem is shown in the figure. The new global 

objective function is a measure of the deviations of the individual 

★ 

objective functions , W and a, from their respective optimum values W and 
a* and is denoted by F (<|>) . 


Multiple objective functions: F.j(<|>) = W 

F 2 (4>) = a 

Constraints, g(<ji): 1-fk/ f k L <0 

V^y-I ? 0 

a - Al < 0 

Global objective function: 



subject to g (<|>)< 0 


FIGURE 19 
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OPTIMIZATION RESULTS FOR RECTANGULAR BLADE USING GLOBAL CRITERIA APPROACH 

(BLADE IN VACUUM) 


Following are the optimization results for the weight-stress optimization 
procedure discussed in the previous chart performed with the blade in 
vacuum. Figure 20 presents results obtained from the single objective 
function compared to those obtained from the multiple objective function 
formulation using the Global Criteria Approach. The results are for the 
rectangular blade with 30 design variables (t^ t 2 and t^ at ten spanwise 

locations) . Case 1 corresponds to the values obtained after optimization 
with blade weight as the single objective function and Case 2 refers to the 
values obtained after optimization with maximum centrifugal stress as the 
single objective function. Case 3 corresponds to the values obtained after 
optimization with multiple objective functions (blade weight and maximum 
centrifugal stress) using the Global Criteria Approach. When only the blade 
weight is minimized, the blade stress increases (Case 1) . On the other hand 
when blade stress is minimized, the blade weight increases (Case 2) . As 
shown using the Global Criteria Approach (Case 3) , when considering both 
stress and blade weight simultaneously, the optimum results fall in between 
those obtained using only single objective functions. Compared to Case 1 
the blade weight is slightly larger but the stress is much lower. Compared 
to Case 2 the blade weight is much lower and the stress is only slightly 
increased. The Global Criteria Approach therefore provides the 'best' 
compromise when two such conflicting objective functions are used. 



Stress 

jV 

H ■ 

***** 


1 

m 



Case 

1 

2 

3 


Case 1 : Objective function 

□ weight 

Case 2: Objective function 
= stress 

Case 3: Objective function 

□ weight & stress 


FIGURE 20 
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CONCLUDING REMARKS 


The paper addresses the problem of structural optimization of helicopter 
rotor blades with integrated dynamic and aerodynamic design considerations. 
Results of recent optimization work on rotor blades for minimum weight with 
constraints on multiple coupled natural flap-lag frequencies, blade 
autorotational inertia and centrifugal stress has been reviewed. A strategy 
has been defined for the ongoing activities in the integrated 
dynamic /aerodynamic optimization of rotor blades. As a first step the 
integrated dynamic/airload optimization problem has been formulated. To 
calculate system sensitivity derivatives necessary for the optimization 
recently developed Global Sensitivity Equations (GSE) are being 
investigated. A need for multiple objective functions for the integrated 
optimization problem has been demonstrated and various techniques for 
solving the multiple objective function optimization are being investigated. 
The method called the 'Global Criteria Approach' has been applied to a test 
problem with the blade in vacuum and the blade weight and the centrifugal 
stress as the multiple objectives. The results indicate that the method is 
quite effective in solving optimization problems with conflicting objective 
functions . 


• Reviewed procedure for dynamic optimization 

with minimum weight objective and frequency, 
autorotational inertia and stress constraints 

• Defined strategy for integrating the above 

with complete aerodynamic optimization 

• Formulated integrated dynamic/airload 

optimization 

• Investigating global sensitivity equations for 

calculating system sensitivity derivatives 

• Described need for multiple objective functions 

• Investigated 'Global Criteria' approach for 

multiple objective optimization 


FIGURE 21 
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